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Abstract 



We calculate bubble-nucleation rates for (2+l)-dimensional scalar theories at high 
^ . temperature. Our approach is based on the notion of a real coarse-grained potential. 

The region of applicability of our method is determined through internal consistency 
criteria. We compare our results with data from lattice simulations. Good agreement 
is observed when the renormalized action of the simulated theory is known. 



1 Introduction 



1.1 The standard formalism 

The standard formalism for the calculation of bubble-nucleation rates during first-order phase tran- 
sitions in field theories at non-zero temperature was introduced in refs. [Q] — [Q] - It consists of an 
implementation of Langer's theory of homogeneous nucleation [|| (for a review see ref. ||) within the 
field-theoretical context. The nucleation rate I gives the probability per unit time and volume to 
nucleate a certain region of the stable phase (the true vacuum) within the metastable phase (the false 
vacuum). Its calculation relies on a semiclassical approximation around a dominant saddle-point, 
which is identified with the critical bubble. This is a static configuration (usually assumed to be 
spherically symmetric) within the metastable phase whose interior consists of the stable phase. It 
has a certain radius that can be determined from the parameters of the underlying theory. Bubbles 
slightly larger than the critical one expand rapidly, thus converting the metastable phase into the 
stable one. The nucleation rate is exponentially suppressed by the action (the free energy rescaled 
by the temperature) of the critical bubble. Possible deformations of the critical bubble generate a 
pre-exponential factor. The leading contribution to this factor has the form of a ratio of fluctuation 
determinants and corresponds to the first-order correction to the semiclassical result in a systematic 
expansion around the saddle point. 

For a (d + l)-dimensional theory of a real scalar field at temperature T, in the limit that thermal 
fluctuations dominate over quantum fluctuations, the bubble-nucleation rate is given by 0] 



Eo(S_\ d/2 
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det^r/^ 2 ]^ 



det^WVo 



exp(-S). (1.1) 



Here T is the free energy of the system for a given configuration of the field <fi. The rescaled free energy 
of the critical bubble is S = F^/T = [F (0b(?")) ~~ ^(0)] /T, where (j)b(r) is the spherically-symmetric 
bubble configuration and = corresponds to the false vacuum. The prime in the fluctuation 
determinant around the bubble denotes that the d zero eigenvalues of the operator [5 2 F /dcj) 2 ]^^, 
corresponding to displacements of the bubble, have been removed. Their contribution generates the 
factor (S/27r) d / 2 and the volume factor that is absorbed in the definition of / (nucleation rate per unit 
volume). The quantity Eq is the square root of the absolute value of the unique negative eigenvalue. 

In field theory, the free energy density T of a system for homogeneous configurations is usually 
identified with the temperature-dependent effective potential. This is evaluated through some per- 
turbative scheme, such as the loop expansion. The profile and the free energy of the critical bubble 
are determined through the potential. This approach, however, faces fundamental difficulties: For 
example, the effective potential, being the Legendre transform of the generating functional for the 
connected Green functions, is a convex function of the field. Consequently, it does not seem to be the 
appropriate quantity for the study of tunnelling. Also, the fluctuation determinants in the expression 
for the nucleation rate have a form completely analogous to the one-loop correction to the potential. 
The question of double-counting the effect of fluctuations (in the potential and the prefactor) must be 
properly addressed. A closely related issue concerns the ultraviolet divergences that are inherent in 
the calculation of the fluctuation determinants in the prefactor. An appropriate regularization scheme 
must be employed in order to control them {jj. Moreover, this scheme must be consistent with the 
one employed for the absorption of the divergences appearing in the calculation of the potential. 

1.2 Coarse-graining 



In refs. UHll! ft was shown that all the above issues can be resolved through the implemention of 
the notion of coarse graining in the formalism. The appropriate quantity for the description of the 
physical system is the effective average action F^ O], which is the generalization in the continuum 



of the blockspin action of Kadanoff [13]. It can be interpreted as a coarse-grained free energy at a 
given scale k. Fluctuations with characteristic momenta q 2 <; k 2 are integrated out and their effect is 
incorporated in T^. In the limit k — ► 0, Tj, becomes equal to the effective action. The k dependence of 
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T/j is described by an exact flow equation [|14| , typical of the Wilson approach to the renormalization 
group [ 15 1 . This flow equation can be translated into evolution equations for the functions appearing 
in a derivative expansion of the action [16, 17]. Usually, one considers only the effective average 
potential Uk and a standard kinetic term, and neglects higher derivative terms in the action. We shall 
employ this approximation in this paper also. The bare theory is defined at some high scale A that can 
be identified with the ultraviolet cutoff. At scales k below the temperature T, a {d + l)-dimensional 
theory at non-zero temperature can be described in terms of an effective (f-dimensional action at zero 
temperature [18]. 

In ref. ||] we considered a (3+l)-dimensional theory of a real scalar field at non-zero temperature, 
defined through its action at a scale ko below the temperature, so that the theory has an effective 
three-dimensional description. The form of the potential Uk results from the bare potential U\ after 
the integration of (quantum and thermal) fluctuations between the scales A and k^. We computed 
the form of the Uk at scales k < k^ by integrating an evolution equation derived from the exact flow 
equation for IV Uk is non-convex for non-zero k, and approaches convexity only in the limit k — > 0. 
The nucleation rate must be computed for k larger than the scale kf at which the functional integral 
in the definition of Uk starts receiving contributions from field configurations that interpolate between 
the two minima. This happens when — k 2 becomes approximately equal to the negative curvature 
at the top of the barrier fl9|| . For k ^ kf the typical length scale of a thick- wall critical bubble is 
<; 1/k. We performed the calculation of the nucleation rate for a range of scales above and near kf, 
for which Uk is non-convex. In our approach the pre-exponential factor is well-defined and finite, as an 
ultraviolet cutoff of order k is implemented in the calculation of the fluctuation determinants, so that 
fluctuations with characteristic momenta q 2 J; k 2 are not included. This is a natural consequence of 
the fact that all fluctuations with typical momenta above k are already incorporated in the form of Uk- 
This modification also resolves naturally the problem of double-counting the effect of the fluctuations. 

We found that the saddle-point configuration has an action Sk with a significant k dependence. 
For strongly first-order phase transitions, the nucleation rate / = ^ exp(— Sk) is dominated by the 
exponential suppression. The main role of the prefactor Ak, which is also k dependent, is to remove 
the scale dependence from the total nucleation rate. Thus, this physical quantity is independent of the 
scale k that we introduced as a calculational tool. The implication of our results is that the critical 
bubble should not be identified just with the saddle point of the semiclassical approximation. It is 
the combination of the saddle point and its possible deformations in the thermal bath (accounted 
for by the fluctuation determinant in the prefactor) that has physical meaning. We also found that, 
for progressively more weakly first-order phase transitions, the difference between Sk and ln(Ak/kf) 
diminishes. This indicates that the effects of fluctuations become more enhanced. At the same 
time, a significant k dependence of the predicted nucleation rate develops. The reason for the above 
deficiency is clear. When the nucleation rate is roughly equal to or smaller than the contribution from 
the prefactor, the effect of the next order in the expansion around the saddle point is important and 
can no longer be neglected. This indicates that there is a limit to the validity of Langer's picture of 
homogeneous nucleation. The region of validity of this picture was investigated in detail in ref. [10|. 
First-order phase transitions in two-scalar models were studied in ref. [11] and the consistency of the 
approach summarized above was reconfirmed. Moreover, the applicability of homogeneous nucleation 
theory to radiatively-induced first-order phase transitions was tested. It was found that the expansion 
around the semiclassical saddle point is not convergent for such phase transitions. This indicates that 
estimates of bubble-nucleation rates for the electroweak phase transition that are based only on the 
saddle-point action may be very misleading. 



1.3 Plan of the paper 

In this paper we present an application of our formalism to (2+l)-dimensional theories at non-zero 
temperature. Our investigation provides a new test of several points of our approach that depend 
strongly on the dimensionality, such as the form of the evolution equation of the potential, the nature 
of the ultraviolet divergences of the fluctuation determinants, and the k dependence of the saddle- 
point action and prefactor. The complementarity between the k dependence of Sk and Ak is a crucial 
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requirement for the nucleation rate / to be k independent. A strong motivation for this study stems 



from the existence of lattice simulations of nucleation for (1+1) and (2+l)-dimensional systems |2(J, 21 
In particular, we shall compare our predictions for the nucleation rate with the lattice results of ref. [21|. 



In the next section we summarize the basic steps of our method and derive the necessary expressions 
for the calculation of the nucleation rate. In section 3 we present sample calculations in two dimensions. 
In section 4 we apply our formalism to theories that have been studied through lattice simulations 
and compare with the lattice results. Our conclusions are presented in section 5. 



2 The calculation of bubble-nucleation rates 



2.1 Evolution equation for the potential 

We consider a model of a real scalar field <f> in 2+1 dimensions. The effective average action Tk((f>) [HU 
is obtained by adding an infrared cutoff term to the bare action, so that contributions from modes with 
characteristic momenta q 2 ^ k 2 are not taken into account. We use the simplest choice of a mass- like 
cutoff term ~ k 2 cj) 2 , for which the perturbative inverse propagator for massless fields is P k {q) ~ q 2 + k 2 . 
This choice makes the calculation of the fluctuation determinants in the pre-exponential factor of the 
nucleation rate technically feasible. Subsequently, the generating functional for the connected Green 
functions is defined, from which the generating functional for the 1PI Green functions can be obtained 
through a Legendre transformation. The presence of the modified propagator in the above definitions 
results in the effective integration of the fluctuations with q 2 ^ k 2 only. Finally, the effective average 
action is obtained by removing the infrared cutoff from the generating functional for the 1PI Green 
functions. 

The effective average action T^. obeys an exact flow equation, which describes its response to 
variations of the infrared cutoff k 113] . This can be turned into evolution equations for the functions 
appearing in a derivative expansion of [16]. In this work we use an approximation which neglects 
higher derivative terms in the action and approximates it by 

(2.1) 



d'x drf +U k {(f>) 

The above action describes the effective two-dimensional theory that results from the dimensional 
reduction of a high-temperature (2+l)-dimensional theory at scales below the temperature. The 
temperature has been absorbed in a redefinition of the fields and their potential, so that these have 
dimensions appropriate for an effective two-dimensional theory. The correspondence between the 
quantities we use and the ones of the (2+l)-dimensional theory is given by 



4>2+l 



T 



U{cf>) 



U 2 +i (02+1, T) 
T 



(2.2) 



In this way, the temperature does not appear explicitly in our expressions. This has the additional 
advantage of permitting the straightforward application of our results to the problem of quantum 
tunnelling in a two-dimensional theory at zero temperature [[|. 

The evolution equation for the potential can be written in the form [14, [H], 10] 



d 

W 2 



i 

8^ 



In 1 + 



k 2 



In 1 + 



k 2 



(2.3) 



For the numerical integration of the above equation we use the algorithms described in ref. [23]. The 
first step of an iterative solution of eq. (|2.3|) gives p2[ 



u [ k\4>) - u£>(o) = UM) - U ko (0) + £ In 



r(i) 
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det [-d 2 + k 2 + U' k \m det[-d 2 + k 2 + U' k [(0)} 
det[-d 2 + k 2 + U' k \ (</>)] det[-d 2 + k 2 + 17» (0)] 



(2.4) 



The scale k\ in the above expression can be chosen arbitrarily anywhere between ko and k, as the 
induced uncertainty for ujp ((f)) corresponds to a higher-order contribution in the iterative procedure. 
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For ki = ho, k — > 0, eq. ( |2.4| ) is a regularized one- loop approximation to the effective potential. Due 
to the ratio of determinants, only momentum modes with k 2 ^ q 2 5s k 2 are effectively included in 
the momentum integrals. The above expression demonstrates the form of ultraviolet regularization 
of fluctuation determinants that is consistent with the cutoff procedure that leads to the evolution 
equation for the potential. An analogous regularization must be used for the fluctuation determinants 
in the expression for the nucleation rate. 



2.2 The nucleation rate 



The calculation of the nucleation rate proceeds in complete analogy to the one in ref. || which we 
outlined in the introduction. We define the theory through the potential at a scale ko below the 
temperature, so that the behaviour is effectively two-dimensional. We then integrate the evolution 
equation down to a scale k <; kf, where we compute the bubble- nucleation rate. In practice, k 2 is 
taken 10% larger than the absolute value of the curvature at the top of the barrier. The potential 
has two minima: the stable (true) one located at eft = cftt, and the unstable (false) one at <j> = (f>f = 0. 
The nucleation rate is exponentially suppressed by the action S k (the free energy rescaled by the 
temperature) of the saddle-point configuration (ftb( r ) that is associated with tunnelling. This is an 
S'0(2)-symmetric solution of the classical equations of motion which interpolates between the local 
maxima of the potential — Uk((f>)- It satisfies the equation 



d 2 <Pb + 1 d4h = ^ 
dr 2 r dr 



with the boundary conditions </>b 
saddle point is given by 



'kVPbh ( 2 - 5 ) 
for r — > oo and d(pb/dr = for r = 0. The action Sk of the 



S k = 2ii 



d<t>b(r) 
dr 



+ U k (Mr))-U k (0) 



rdr = S{ + S v k . 



(2.6) 



The profile of the saddle point can be easily computed with the "shooting" method [ 24 ] . A consistency 
check for our solution is provided by the fact that S k ! = for two-dimensional theories. 

The bubble-nucleation rate is determined through eq. ( |1.1| ) in terms of the potential Uk(4>)- The 
explicit temperature dependence is absorbed in the definition of effective two-dimensional parameters 
according to eq. ( |2.2D . An appropriate regularization is implemented in order to control the ultraviolet 
divergence of the prefactor. The type of regularization is dictated by the one-loop effective potential, 
given by eq. ([2.4]) in our scheme. This equation indicates that fluctuation determinants computed 
within the low-energy theory must be replaced by appropriate ratios of determinants. We emphasize 
that the matching of the regularization scheme for the prefactor with the cutoff used for the derivation 
of the evolution equation for the potential ( p.3|) is crucial for the consistency of our method. Finally, 
the nucleation rate in two dimensions is given by 



A k exp(-Sfe) 



where 



2vr 2vr 



det' [-d 2 + U'k\4>b{r))] det [-d 2 + k 2 + U' k \0)} 



det [-d 2 + k 2 + UjKMr))] det [-8 2 + C^'(0)] 
The differential operators that appear in eq. ( |2.7| ) have the general form 

W Ka = -d 2 + ml + aW k {r) 



-1/2 



where 



m 2 



U' k \0) + Kk 2 , 
W k (r) = Unfair)) ~ Uk(0), 



(2.7) 



(2.8a) 

(2.8b) 
(2.8c) 
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with k, a = or 1. As the W Ka operators are S'0(2)-symmetric, it is convenient to use polar coordinates 
and express their eigenfunctions as ip n (r,{p) = e mv u n (r)/^/r. This leads to 



detW K 



detHW 



dr 2 



+ 



n 



ml + aW k (r). 



(2.9) 



The computation of such determinants is made possible by a theorem [25] that relates ratios of 
determinants to solutions of ordinary differential equations. In particular, we have 



detW 



9n 



nul 



oo 



detW ntt0 VnKo{r 
where y n na(r) is the solution of the differential equation 

d 2 n 2 



oo 



(2.10) 



dr 2 



+ 



+ m 2 K + aW k {r) 



Unnair) = 0, 



(2.11) 



with the behaviour y nKa (r) oc r'™' + 2 for r — > 0. For example, y nK o are proportional to modified 
Bessel functions: y nK o oc ^fr I\ n \{m K r). Equations such as ( |2,11 ) can be solved numerically with 
Mathematica [p6f] . Since opposite values of n lead to identical determinants, the final expression for 
the nucleation rate can be written as 



co 



i(l) exp( - Sfc) n o cn ' 

1/2 



1 9oo I 



Cl 



9u_ 



9nl 
9n0 



(2.12) 



The calculation of c\ is slightly complicated because of the necessity to eliminate the zero eigenvalue 
in g' 10 . (The two zero eigenvalues of the operator — d 2 + U' k \(f)\>{r)) are included in the equal factors 
<7_io and gio). Also the (unique) negative eigenvalue — E 2 must be computed for the determination of 
Co- How these steps are achieved is described in ref. ||. For sufficiently large re, one can compute c n 
analytically using first-order perturbation theory in Wk(r) 0, ||. We find 



9,, 



1 + 



2re 



r Wk(r) dr, 



1 



4re 3 



fc / r 6 W k (r)dr. 



(2.13) 



These expressions are very useful for the evaluation of the prefactor, as only c n for small values of n 
need to be computed numerically. 



3 A sample computation 

We are interested in potentials that have the approximate form 



(3.1) 



The explicit temperature dependence has been absorbed in 7, h and <f> according to eqs. (| 
a shift <b — ► c 



. Through 

c the cubic term in eq. (3.1) can be eliminated in favour of a term linear in (p. The 



resulting potential describes a statistical system of the Ising universality class in the presence of an 



external magnetic field. We assume that the potential of eq. (3.1) describes the theory at some 
initial scale ko not far below the temperature, similarly to refs. ||]— [11]. Its form at lower scales can 
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Figure 1: The steps in the computation of the nucleation rate: (a) Potentials Vk((ft); (b) Saddle points 
(j)b{r); (c) Wk(r), given by eq. ( 2.8 c\) ; (d) Results for the saddle-point action Sk (diamonds), prefactor 
In (A/c/jn 3 ) (stars) and nucleation rate — In (//to 3 ) (squares). The first row corresponds to a model 
with a potential Uk ((ft) given by eq. (3.1) with "f/m 2 = —2, h/m 2 = 1/3, while the second one to a 
model with 7/m 2 = —4.5, h/m 2 = 2. All dimensionful quantities are given in units of m. 



be determined by integrating the evolution equation (^^) numerically, or by using the approximate 
solution of eq. (^J). 

The various steps in our calculation are summarized in fig. 1, for two models described by a 



potential Uk {(ft) given by eq. (3J). In the following we express all dimensionful quantities in terms 
of the arbitrary mass scale m. In the first row we present results for a theory with 7/771 2 = —2, 
h/m 2 = 1/3. In (a) we present the evolution of the potential Uk((ft) as the scale k is lowered. We 
always shift the metastable vacuum to (ft = 0. The solid line corresponds to k$/m = 2, while the line 
with longest dashes (that has the smallest barrier height) corresponds to kf/m = 1.1. At the scale kf 
the negative curvature at the top of the barrier is slightly larger than —k 2 . This is the point in the 
evolution of the potential where configurations that interpolate between the minima start becoming 
relevant in the functional integral that defines the coarse-grained potential [19|. For this reason, we 
stop the evolution at this point. We observe that, for low k, the absolute minimum of the potential 
settles at a non-zero value of (ft. A significant barrier separates it from the metastable minimum at 
= 0. The profile of the saddle point <^>b( r ) is plotted in (b) in units of m for the same sequence 
of scales. We observe a variation of the value of the field (ft in the center of the critical bubble for 
different k. This is reflected in the form of the quantity Wf-(r), defined in eq. Q2.8c| ), which we plot 
in (c). 

Our results for the nucleation rate are presented in (d). On the horizontal axis we give the values 
of ln(fc/m). The dark diamonds correspond to the values of the action S% of the saddle point at 
the scale k. The stars indicate the values of ln(A^/m 3 ). We observe a logarithmic k dependence for 
both these quantities. The value of is expected to decrease for smaller k, because k acts as the 
effective ultraviolet cutoff in the calculation of the fluctuation determinants in A^. For smaller k, 



6 



fewer fluctuations with wavelengths above an increasing length scale ~ 1/k contribute explicitly to 
the fluctuation determinants. The logarithmic dependence on k is the reflection of the logarithmic 
ultraviolet divergence of the unregularized prefactor in two dimensions. The dark squares in (d) give 
our results for — ln(I/m 3 ) = Sk — ln(Ak/m 3 ). The logarithmic k dependence largely cancels between 
<S/% and ln(y4fc/m 3 ), so that ln(7/m 3 ) is almost constant. The small residual dependence on k can be 
used to estimate the contribution of the next order in the expansion around the saddle point. This 
contribution is expected to be smaller than ln(Afc/m 3 ). This behaviour confirms that the nucleation 
rate should be independent of the scale k that we introduced as a calculational tool. 

In the second row we present the calculation of the nucleation rate for a model with a larger 
coupling h/m? = 2, and 7/rra 2 = —4.5, ko/m = 2, kf/m = 1.1. We observe a more pronounced k 
dependence of the potential, saddle-point profile and function Wk{r). The most important aspect of 
the comparison of the two models concerns the relative values of Sk and ln(Ak/m 3 ). For the first 
model, the contribution of the prefactor to the nucleation rate is much smaller than that of the action 
of the saddle point. The main role of the prefactor is to remove the logarithmic k dependence from 
I/m 3 . For the second model, Sk and ln(Ak/m 3 ) are comparable. This indicates that the effects of 
fluctuations are enhanced. Moreover, the prefactor fails to cancel the k dependence of the saddle-point 
action. The reason is that the next order in the expansion around the saddle point is important and 
can no longer be neglected. This establishes the limit of validity of homogeneous nucleation theory || 
for two-dimensional systems, in agreement with the studies of refs. 0]-[ll| in three dimensions. 



4 Comparison with lattice studies 
4.1 Matching the lattice theory 

A main objective of this work is the comparison of our results with the lattice study of ref. pi]] . This 
requires a precise definition of the form of the potential. We must make sure that we consider a 
theory identical to the one simulated on the lattice. As we cannot match exactly the bare parameters 
of the lattice action, it is more convenient to guarantee that the low energy renormalized theory is 
the same in our model and ref. [21|. In the latter work, through a redefinition of the field and the 



distance, the one-loop renormalized action (free energy rescaled with respect to the temperature) of 
the (2+l)-dimensional theory at high temperature is expressed as 



with 



(1 " 4^) W + v$) + — [v"(4>) - V"(4>)\n (v"(4>))} }, (4.2) 

V(4>) = \tf - \^ + ^A0 4 . (4.3) 

For A < the potential becomes unbounded from below, while for A = 1/3 it has two equivalent 
minima. For this reason, we consider A values in the interval (0,1/3). Higher- loop corrections are 
expected to be proportional to larger powers of 9/8tt. This implies that the theory that is simulated 
on the lattice corresponds to a renormalized action in the continuum given by eqs. ( |4.2| ), fl4.3| ) only 
if 6 is not much larger than 1. In the opposite case, the renormalized action cannot be determined 
perturbatively. 

In our approach the theory is defined at some initial scale ko and the integration of evolution 
equations such as eq. (|2.3| ) generates the low-energy structure. For k = the effective average action 



of eq. (2A) must be identified with the action of eq. ( ^4.2[) . We neglect wavefunction renormalization 
effects, which can be seen from eq. ( |4.2| ) to be a good approximation for values of 9 not much larger 
than 1. The form of the potential at some non-zero scale k can be inferred by demanding that the 



iterative solution of eq. (2.4) reproduces the potential of eq. ( [4.2[ ) for k = 0. In particular, by choosing 



k = k\ = in eq. (|2.4|) and renaming ko as k we find 



Um - VW + ^V"{<t>) ~l(k 2 + V»W) In " + J 2 //(0) , (4.4) 
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where V((j>) is given by eq. (3.1) with 



4, =->/*, —2=\^- (4-5) 
m z m z 3 

In eq. ( ^4.4| ) we have neglected terms ~ (87r) -2 and terms ~ (87r) _1 in the arguments of the logarithms. 
It is clear from the above that the dimensionless coupling that controls the validity of the perturbative 
expansion is 9X/3. For this reason, perturbation theory is expected to break down for 9 J; 3/ A. 

In summary, eq. (fO|) is an approximate solution of the evolution equation (|2.3|) (at the first 
level of an iterative procedure) which is consistent with eq. (^]^) that determines the renormalized 
action of ref. The matching between the renormalized and the lattice actions is accurate only at 
the one- loop level (the region of validity of eq. Q4.2| )). The approximate solution of eq. (|Q|), which 
results at the first level of the iterative procedure, has a similar region of applicability. Morever, this 
approximation is not valid for values of k that render negative the argument of the logarithm in the 
right-hand side of eq. (|4.4|), i.e. for k 2 < max {— V" (</>)}• This implies that eq. (4.2) is trustable only 



in the convex regions of the potential (for which V"{4>) > 0) and should not be expected to lead 



to reliable predictions for the nucleation rate. This is in agreement with the conclusion of ref. [21 



that the determination of the bubble- nucleation rate through the real part of eq. (4.2) does not lead 
to consistency with the lattice results. On the other hand, eq. ( (4.2j ) is perfectly valid in the convex 
regions of the potential, where it can be matched with eq. (|4.4[) for k = 0. In the following, for our 



predictions of the bubble-nucleation rates and the comparison with the lattice results, we rely on the 



approximate solution of eq. (|4.4| ), instead of integrating eq. numerically. The numerical solution, 
which is more accurate than eq. ( ^4.4[) , does not offer increased precision in our comparison with the 
lattice data, as the determination of the renormalized theory is valid only at the one-loop levelQ. 

4.2 Comparison with the lattice results 

In fig. 2 we present a comparison of results obtained through our method with the lattice results of 
fig. 1 of ref. [^] . For each of several values of A we vary the parameter 9 and determine the couplings 
7, h according to eqs. Q4,5| ). The coarse-grained potential is then given by eq. ( fj.4j ) for k > kf. The 
diamonds denote the saddle-point action Sk- For every choice of A, 9 we determine Sk at two scales: 
l-2kf and 2 k f. The light-grey region between the corresponding points gives an indication of the k 
dependence Sk- The bubble-nucleation rate — ln(//m 3 ) is denoted by dark squares. The dark-grey 
region between the values obtained at 1.2 kj and 2kf gives a good check of the convergence of the 
expansion around the saddle-point. If this region is thin, the prefactor is in general small and cancels 
the k dependence of the action. The dark circles denote the results for the nucleation rate from 
the lattice study of ref. J21[| . The dashed straight lines correspond to the action of the saddle-point 
computed from the 'tree-level' potential of eq. Q3,l|). 



For A = the potential of eq. (3.1) is unbounded from below and the procedure we outlined in the 
previous paragraphs for matching the theory simulated on the lattice may be problematic. However, 
one may consider these results as applying to the limit A — > 0, so that no conceptual problems arise. 
The values of — ln(//m 3 ) computed at 1.2 kf and 2 kf are equal to a very good approximation, which 
confirms the convergence of the expansion around the saddle-point and the reliability of the calculation. 
The k dependence of the saddle-point action is cancelled by the prefactor, so that the total nucleation 
rate is k independent. Moreover, the prefactor is always significantly smaller than the saddle-point 



action. The circles indicate the results of the lattice simulations of ref. |21]. The agreement with the 
lattice predictions is good. More specifically, it is clear that the contribution of the prefactor is crucial 
for the correct determination of the total bubble-nucleation rate. Similar conclusions can be drawn 
for A = 0.1 and A = 0.2. 

For larger values of A the lattice simulations have been performed only for 9 significantly larger 
than 1. For smaller 9, nucleation events become too rare to be observable on the lattice. As we 



1 We also perform checks of the corrections arising from integrating the full evolution equation (2.3). These corrections 



are very small in the parameter region for which the expansion around the saddle point is convergent. 
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Figure 2: Comparison of our method with lattice studies: Diamonds denote the saddle-point action 
Sk and squares the bubble-nucleation rate — In (I/m 3 ) for k = 1.2 kt and 2kf. Dark circles denote 
the results for the nucleation rate from the lattice study of ref. Finally, the dashed straight lines 

correspond to the action of the saddle-point computed from the potential of eq. ( \3. \ ). 



discussed earlier, the matching between the lattice and the renormalized actions becomes imprecise 
for large 9. This indicates that we should expect deviations of our results from the lattice ones, as 
the theory of eq. (|4.2| ) may be different than the simulated one. These deviations start becoming 
apparent for the value A = 0.25, for which the lattice simulations were performed with 8 ~ 10-20. For 
A > 0.3 the lattice results are in a region in which the internal consistency criteria of our method for 
the reliability of the expansion around the saddle-point are not satisfied any more. 

The consistency of our calculation is achieved for 1/0 ^ 0.12 even for A = 0.32. However, the 
breakdown of the expansion around the saddle point is apparent for 1/9 0.12. The k dependence 
of the predicted bubble-nucleation rate is strong^. The prefactor becomes comparable to the saddle- 
point action and the higher-order corrections are expected to be large. The k dependence of S k is 
very large. For this reason we have not given values of Sk in this case. This behaviour indicates that 
the field fluctuations become significant. Typically, these fluctuations enhance the total rate and for 
1/9 ^ 0.08 can even compensate the exponential suppression. Several similar examples were given in 
refs. g-Q. 

The reason for this behaviour can be traced to the form of the differential operators in the prefactor 
(see eqs. eqs. (|2.7| )-( ^8c]) ). This prefactor, before regularization, involves the ratio det'(— d 2 -\-U' k '(0) + 
W k (r))/ det(-d 2 + U' k \0)), with W k (r) = U'^((j) b (r)) - U' k \0). The function W k (r) always has a 

2 The additional squares for 1/6 = 0.1 correspond to results from the numerical integration of eq. (^.3|). 
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minimum away from r = (see figs, lc and lg), where it takes negative values. As a result the lowest 
eigenvalues of the operator det'(— d 2 + U£(0) + Wfc(r)) are smaller than those of det(— d 2 + U^(0)). 
The elimination of the very large eigenvalues from the determinants through regularization does not 
affect this fact and the prefactor is always larger than 1. Moreover, in cases such as those depicted 
in fig. 2f it becomes exponentially large because of the proliferation of low eigenvalues in det'(— d 2 + 
U'f!(0) + Wk(r)). In physical terms, this implies the existence of a large class of field configurations of 
free energy comparable to that of the saddle-point. Despite the fact that they are not saddle points of 
the free energy (they are rather deformations of such a point) and are, therefore, unstable, they result 
in a dramatic increase of the nucleation rate. This picture is similar to that of "subcritical bubbles" 
of ref. [27]. In ref. [28] the nucleation rate was computed by first calculating a corrected potential 
that incorporates the effect of such non-perturbative configurations. The pre-exponential factor must 
be assumed to be of order 1 in this approach, as the effect of most deformations of the critical bubble 
has already been taken into account in the potential. In our approach the non-perturbative effects are 
incorporated through the prefactor. Both methods lead to similar conclusions for the enhancement of 
the total nucleation rate. 

A final comment concerns the 9 dependence of our results for the nucleation rate in fig. 2. For 
a given value of A, we observe a linear dependence of — ln(I/m 3 ) on 1/9 in the regions where the 
calculation is consistent. This is explained by the fact that the dominant 9 dependence of the potential 
arises from the first term in the right-hand side of eq. (|4.4j ). This is in agreement with the findings of 
ref. pjjl . I n the latter work, this linear dependence was not observed for A = 0.3 and 0.32. The reason 
is that the range of 9 used in the simulations was so large that the first term in the right-hand side 
of eq. fl4.4Q ceased to be dominant. Our results for small 9 display the expected behaviour even for 
A = 0.32. 



5 Conclusions 

In this paper we applied our approach for the calculation of bubble-nucleation rates to (2+l)-dimen- 
sional theories at non-zero temperature. We studied these theories at coarse-graining scales k below 
the temperature, where they display an effective two-dimensional behaviour. This provided the oppor- 
tunity to check several points of our approach that depend on the dimensionality of the system. For 
example, the evolution equation for the coarse-grained potential has a different form than the one in 
the effective three-dimensional systems we studied in the past. Also, the pre-exponential factor in the 
bubble-nucleation rate has a logarithmic ultraviolet divergence before regularization, instead of the 
leading linear divergence in three dimensions. This divergence is reflected in the leading logarithmic 
dependence of the regularized prefactor on the scale k that acts as an ultraviolet cutoff (fig. 1). More 
crucially, the complementarity between the k dependence of the prefactor and the saddle-point action 
was observed again. This is a crucial point that guarantees that a physical quantity, such as the 
bubble-nucleation rate, is independent of the scale k that we introduced as a calculational tool. The 
validity of Langer's theory of homogeneous nucleation was confirmed, as long as the prefactor gave a 
contribution to the nucleation rate smaller than the leading exponential suppression by the action of 
the saddle point. 

Another important aspect of this work concerns the comparison of our results with data from 
lattice simulations. This constitutes a stringent quantitative test of our method. We found good 
agreement between our results and the lattice data when the renormalized action for the theory that 
is simulated on the lattice is known (fig. 2). In these cases, we first match the renormalized action and 
then compute the nucleation rate. For part of the range of the lattice parameters this is not possible, 
because the renormalized action cannot be obtained from the lattice action perturbatively. However, 
the internal consistency criteria of our method provide a test of the reliability of our results in all 
cases. 
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